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Abstract 

We study the 3-state hexagonal-lattice Potts antiferromagnet by a Monte 
Carlo simulation using the Wang-Swendsen-Kotecky cluster algorithm. We 
study the staggered susceptibility and the correlation length, and we confirm 
that this model is disordered at all temperatures T > 0. We also measure the 
ground-state entropy density. 



1 Introduction 



Antiferromagnetic Potts models JT], |], H are much less well understood than their 
ferromagnetic counterparts. One reason for this is that the behavior depends strongly 
on the microscopic lattice structure, in contrast to the universality typically enjoyed 
by ferromagnets. As a result, many basic questions have to be investigated case- by- 
case: Is there a phase transition at finite temperature, and if so, of what order? What 
is the nature of the low-temperature phase? If there is a critical point, what are the 
critical exponents and the universality classes? Can these exponents be understood 
(for two-dimensional models) in terms of conformal field theory? 

One thing is known rigorously @, [§]: for q large enough (how large depends on the 
lattice in question), the antiferromagnetic g-state Potts model has a unique infinite- 
volume Gibbs measure and exponential decay of correlations at all temperatures, 
including zero temperature: the system is disordered as a result of the large ground- 
state entropy. However, for smaller values of q, phase transitions can and do occur. 
One expects that for each lattice C there will be a value q c (C) such that 

(a) For q > q c (C) the model has exponential decay of correlations uniformly at all 
temperatures, including zero temperature. 

(b) For q = q c (C) the model has a critical point at zero temperature. 

(c) For q < q c {£) any behavior is possible. Often (though not always) the model 
has a phase transition at nonzero temperature, which may be of either first or 
second order. 

The problem, for each lattice, is to find q c (C) and to determine the precise behavior 
for each q < q c (C). 

The g-state Potts model on a lattice C is defined by the Hamiltonian 

n = -J £ a- 1 ) 

where the sum J2(xy > runs over & U possible nearest-neighbor pairs of lattice sites (each 
pair counted once), and each spin takes values o~$ G {1, . . . , q}. The antiferromagnetic 
case corresponds to J = —f3 < 0. The Potts antiferromagnet on the hexagonal 
(honeycomb) lattice is the simplest case. First, we know that the q = 2 model has 
a critical point at e J = 2 — V3 because the corresponding ferromagnetic model is 
critical at e J = 2 + y/3, and the hexagonal lattice is bipartite. Second, the results 
of || show that the hexagonal-lattice antiferromagnet exhibits exponential decay of 
correlations at zero temperature for q > 4. Thus, q c {hc) should lie somewhere strictly 
between q = 2 and q — 4. The behavior at q = 3 has not yet been analyzed rigorously. 

On the other hand, the g-state hexagonal-lattice Potts model has been solved 
(e.g., the free and internal energies are exactly known) on two special curves in the 
(J, q) plane (see Fig. [1]): 

(e J -l) 3 -3g(e J -l) = q 2 (1.2) 
e J = l-q for < q < 4 (1.3) 



2 



These curves are the duality images |J of the corresponding curves for the triangular- 
lattice Potts antiferromagnet obtained by Baxter and collaborators |7], || ||. Curve 
( |1.2| ) has three branches in the region q > 0; the uppermost branch (with < q < oo 
and e J > 1) corresponds to the ferromagnetic critical line; the middle branch (with 
< q < 4 and — 1 < e J < 1) contains the q = 2 antiferromagnetic critical point 
(e J = 2 — a/3) and crosses the zero-temperature line e J = at q = (3 + \/5)/2 « 2.618; 
the lower-most branch crosses the zero-temperature line at q = (3 — Vo)/2 ~ 0.382.Q 
The second curve ( |1.3|) is physical (e J > 0) only for < q < 1. 

The behavior of the middle branch suggest that it may be the antiferromagnetic 
critical curve. If this is the case, there would be a zero-temperature critical point at 

q = q c {hc) = (1.4) 

(if this assertion has any meaning^) and the model would be disordered at all tem- 
peratures for q > q c (hc). This conclusion is in agreement with the exact result of 
Ref. ||. Further interesting speculations can be found in Ref. |TD . 



Shrock and Tsai |15[ used this theoretical argument to rule out the existence 
of a critical point in the 3-state antiferromagnet. They also pointed out that this 
argument could fail if there were a first-order phase transition (with finite correlation 
length). To test these ideas, they performed Monte Carlo simulations over a large 
range of temperatures (0 < (3 < 5), and studied the (static and dynamic) behavior of 
the energy density. They found no signal of hysteresis, and no critical slowing-down 
at large (3 for the Metropolis |]TB| and Swendsen-Wang cluster |T7|, |TS| algorithms. 



They finally concluded that the model is disordered at all temperatures T > 0. 

However, the numerical results of Ref. JTSJ (namely, smoothness of the energy 
density and boundedness of its integrated autocorrelation time for Metropolis and 
Swendsen-Wang algorithms) do not constitute a strong evidence supporting the ab- 
sence of criticality of this model. A closely related model (e.g., the 3-state square- 
lattice Potts antiferromagnet || plj) makes an excellent counter-example: it has a 



critical point at T = 0, but the energy density is smooth over the entire temperature 
range, and the autocorrelation times for the Wang-Swendsen-Kotecky (WSK) cluster 
algorithm JZUj, ^T| are uniformly bounded r int < 7.6 |22|, ^3|, On the other hand, 



the absence of critical slowing-down at large f3 for the (local) Metropolis algorithm 
only implies that the specific heat Ch is bounded. This comes from the rigorous 



The middle branch is missing in Ref. @, p. 673, Figure 8. 

2 The Potts models for non-integer q can be given a rigorous meaning via the mapping onto the 
Fortuin-Kasteleyn random-cluster model |ll|, [l2| The trouble is that in the antiferromagnetic 
case ( J < 0) this latter model has negative weights, and so cannot be given a standard probabilistic 
interpretation. In particular, the existence of a good infinite- volume limit is problematical; the limit 
could depend strongly on the subsequence of lattice sizes and on the boundary conditions. The 
same is true of the "anti-Fortuin-Kasteleyn" representation, in which the coefficients are products 
of chromatic polynomials of clusters: again the weights can be negative for non-integer q, and the 
existence of the infinite- volume limit is problematical. Likewise, the ice-model representation |g[ Q 
has in general complex weights for < q < 4, even in the ferromagnetic case. 
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bound 

r- m t, £ ZVC H , (1.5) 

where V is the volume. Thus, a direct test of non-criticality for the 3-state hexagonal- 
lattice Potts antiferromagnet is still lacking. This test can be achieved, for instance, 
by measuring the staggered susceptibility and the correlation length. 

Finally, the above-mentioned argument, which identifies exact solubility with criti- 
cality, though plausible, is not necessarily valid. For example, in the triangular-lattice 
case there are two curves where the model can be solved, and the antiferromagnetic 
critical point of the q = 2 model (namely, e J = 0) happens to lie on both of these 
curves. Moreover, the antiferromagnetic critical point of the q = 4 model, which is 
believed || |26], |27|, |28|f] to be at e J = 0, lies on one (though not the other) of these 



curves. Nevertheless, neither of these two curves can properly be identified with the 
antiferromagnetic critical curve of this model, as this identification would predict an 
incorrect scenario for q = 3 (see Ref. || for details). 

In this note, we report the results of performing a direct test of non-criticality and 
we show with no ambiguities that the 3-state Potts antiferromagnet on the hexagonal 
lattice is always disordered. In Section |] we describe the method we have used to 
simulate the system and the operators we have measured. In Section ^ we display 
the results for the energy, specific heat, staggered susceptibility and second-moment 
correlation length. All these quantities (as well as the corresponding integrated au- 
tocorrelation times) are bounded uniformly in the temperature and the lattice size. 
We conclude that the 3-state hexagonal-lattice Potts antiferromagnet is disordered 
at all temperatures T > 0. As a by-product of our calculation, we compute the 
zero-temperature entropy density. 



2 Numerical Simulations 

We have performed Monte Carlo simulations of the 3-state hexagonal-lattice Potts 
antiferromagnet at temperatures ranging from T = oo to T = 0. More precisely, we 
have simulated this model from (3 = to (3 = 9 in intervals of 0.05, and also exactly 
at (3 = oo. We have made our simulations using the WSK cluster algorithm. This 
algorithm is suitable to simulate this model at zero temperature because our lattices 
are bipartite p9| , 24|. At (3 = we have started the simulations with a random 



configuration; at finite (3 > 0, we started from the last configuration generated at the 
closest smallest (3; and at (3 = oo we started from an ordered configuration (spins in 
one sublattice all equal to 1, and spins on the other sublattice all equal to 2). 



The hexagonal lattice is not a Bravais lattice |p0|| , as not all points are equivalent. 
Rather, it is the union of two sublattices, the even and the odd, each of which is 
isomorphic to a triangular lattice (whose lattice spacing is larger by a factor of y/3). 

3 The authors of Ref. [^7| studied the 3-coloring model on the hexagonal lattice, which is equiv- 
alent to the 3-state Potts antiferromagnet on the Kagome lattice at zero temperature. This latter 
model can be exactly mapped onto the 4-state Potts antiferromagnet on the triangular lattice at 
zero temperature fl28fl. 
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The hexagonal lattice can thus be viewed as an underlying triangular lattice (which 
is Bravais) with a two-point basis. To be more precise, consider a finite hexagonal 
lattice Ti with periodic boundary conditions. Then the even sublattice of TC is a 
triangular lattice T. Conversely, given an L x L triangular lattice T with periodic 
boundary conditions, we can construct a hexagonal lattice TC with Vh c = 2L 2 points 
by taking the sites in T together with the centers of the down-pointing elementary 
triangles of T. Thus, a generic point x of the hexagonal lattice can be written as 







x = x'xffi + x' 2 f)2 + eff = x' + erf (2.1a) 
x[,x' 2 = 1,...,L (2.1b) 

where x' lives on the triangular lattice spanned by the (unit) vectors 

4 = (1,0) (2.2a) 

and L is the linear size of the triangular sublattices.Q In this paper we have considered 
lattices ranging from L = 3 to L = 48. The variable e = 0, 1 can be interpreted as 
the "parity" of the corresponding lattice site: if e = (resp. e = 1) the site x belongs 
to the even (resp. odd) sublattice. The vector 

V = ^(0,1) = 7^i + 2t7 2 ) (2.3) 

is the so-called basis vector joining the two points of the basis. The pair (x',e) 
determines uniquely a point on the hexagonal lattice, and conversely, given a point x 
of the hexagonal lattice we can uniquely obtain the pair (x',e) associated to it. 
We have measured three basic observables. The simplest one is the energy 

£ = E • (2-4) 

The staggered magnetization can be written easily if we represent the spin at site x 
by a vector in the plane 

/ 2tt 2tt \ . . 

a s = I cos -g-eraf, sm —ag J . (2.5) 

In this case, the staggering assigns a phase e l7re = (— l) e depending solely on which 
sublattice the spin is located. The square of the staggered magnetization can be 

4 We have chosen the lattice spacing of the triangular sublattice to be or = 1; the lattice spacing 
of the corresponding hexagonal lattice is therefore a-n = \j \f?>. 

5 This choice is motivated by what happens in the q = 2 case: at T — all the spins on the even 
sublattice take one value (say, 1), and all the spins on the other sublattice take the other value (say, 
2). The natural staggering corresponds to assign a phase e l7r to all the spins on the odd sublattice. 
We can generalize this to the q — 3 case by assigning a general phase e 1 ^ to all the spins on the odd 



sublattice. Then, the contribution of the six smallest momenta ( p.lO[ ) to the observable ( 2.11 ) is the 



same in average only when <j) — (uniform magnetization) , and <f> — n (staggered magnetization) . 
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written as 



stagg 



XX _1 ) £<? z'+er? 



o 3 

-Y 



E( "*-) 'W+eif.a 



(2.6) 



This is a "zero-momentum" observable. In order to estimate the second-moment 
correlation length we have to define the corresponding "smallest-nonzero-momentum" 
observable. The translational invariance of the hexagonal lattice 7i is given by the 
underlying triangular (Bravais) lattice T . Thus, the allowed momenta are those of a 



triangular lattice of size L x L with periodic boundary conditions |30 

2tt 



k -- 
m 1 ,m 2 = 

The momenta are given in the basis 

Pi z 

P2 

defined by the relations 



-mxpi + —m 2 p2 
1,...,L 



°=(o,i), 



ffi ■ Pj = Sij . 



The smallest nonzero momenta are 



2vr^ 2tt^ 2tt ^ 



(2.7a) 
(2.7b) 



(2.8a) 
(2.8b) 

(2.9) 
(2.10) 



all having \k\ = 47r/(v / 3£)- Thus, the smallest-nonzero- momentum observable asso- 
ciated to Q2.6Q is 



3 1 



stagg 



X 



2 6 



6 3 

EE 

n=l a=l 



£(-i)v*^' + ^w 



(2.11) 



The contributions of the wavevectors k n and — k n are exactly the same; thus Eq. fl2.11|) 
can be simplified: 



stagg 



3 1 
- x - 

2 3 



E 

a=l 



+ 



+ 



E 

x' 

E 



E 



2Trix[/L 



3 27ri(o:i-4)/I, 



-27TJ 



<T s l +il ,a 



(2.12) 
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From these observables we have computed the following expectation values: the 
energy density (per spin) 

E=±-(S), (2.13) 

Vhc 

the specific heat 

C H = ^~ var(£) = ±- {(£ 2 ) - (£)*} , (2.14) 

Vhc Vhc 

the staggered susceptibility 

Xstagg = ^(^L ffl ) , (2-15) 
'he 

and the second-moment correlation length 

t _ {Xstagg/ F st ggg — 1) ^ - 

^~ 2vrsin(7r/L) ' ^ ibj 



where F stagg is given by 

Fstagg = Tl—{3~stagg) ■ (2-17) 
Vhc 

For each observable O discussed above we have measured its integrated auto- 
correlation time T intj e> using a self-consistent truncation window of width 6ri nt P2| , 
Appendix C]. 

We have made 5 x 10 5 (resp 3.5 x 10 5 ) iterations for L = 3, 6 (resp. L > 9). 
We have discarded a 10% of them to allow the system to reach equilibrium. The 
autocorrelation times for all the observables were uniformly bounded in /3 and L: 

Tint $> 4 • (2.18) 

For L > 9 there is a sharper bound: r int ^3.0 This means that we have made ~ 
10 5 T int measurements, and we have discarded ~ 10 4 r- int iterations as a (conservative) 
prevention against the existence of any slower mode not considered here. 

We have made our simulations on two Pentium machines at 166 MHz. Each WSK 
iteration took approximately 3 x (2L 2 ) /isec; the total CPU used was approximately 
13 days. 



3 Data Analysis 

In this section we perform all fits using the standard weighted least-squares method. 
As a precaution against corrections to scaling, we impose a lower cutoff L > L min 
on the data points admitted in the fit, and we study systematically the effects of 
varying L m i n on both the estimated parameters and the \ 2 ■ m general, our preferred 

6 The fact that r m t is larger for smaller lattices can be understood because the correlation length 
satisfies £ < 3 for all T and L. For small L, the ratio £/L is not much smaller than 1; however, for 
large L, \/L < 1. 
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fit corresponds to the smallest L min for which the goodness of fit is reasonable (e.g., 
the confidence levelQ is > 10-20%) and for which subsequent increases in L min do not 
cause the x 2 to drop vastly more than one unit per degree of freedom. 

Let us first consider the second-moment correlation length £ = £((3,L) (See Fig- 
ure 0). We see that this observable is, for fixed L, a non-decreasing function of (3 
which asymptotically tends to the value at (3 = oo [i.e., = oo,L)). At fixed /?, 
the function •) is also non-decreasing. For L > 12, the values of £(/?, L) collapses 
well onto a single curve. Furthermore, L) < 3.2 uniformly in (3 and L. Thus, the 
correlation length stays bounded even at T = 0; this observation implies that there 
is no critical point for this model at any temperature T > 0. 

If we fit the values of £(oo, L) to a constant [=£(oo, oo)] we obtain a good fit for 
T — 1 9 

J-'mtn - Lz '; 

f(oo, oo) = 3.0828 ± 0.0098 (3.1) 

with x 2 = 161 (2 DF, level = 45%). Thus, £(/3,L) < f(oo, oo) uniformly in /3 
and L. The numerical value of £(oo,oo) is consistent with the observation that the 
thermodynamic limit is attained in practice (i.e. £(/?, L) <C L) for L > 12. Therefore, 
we do not have to consider larger lattices (i.e., L > 48). 

This scenario also applies to the staggered susceptibility (see Figure |). This 
observable \Xstagg(/3, L)] is also a non-decreasing function of (3 at fixed L, and of L 
at fixed (3. Moreover, for L > 12 the measurements collapse well onto a single curve, 
which is uniformly bounded: Xstagg(f3, L) ^ 20. There is no signal of second-order or 
first-order phase transition at any temperature. 

If we fit the value of the staggered susceptibility at zero temperature Xsta gg (P = 
oo, L) to a constant [= Xstaggi^o, oo)], we obtain a good fit for L min = 24 

Xstaggi^, oo) = 19.070 ± 0.043 (3.2) 

with x 2 = 0.48 (1 DF, level = 49%). 

The energy and specific heat are both non-increasing curves which tend asymp- 
totically to zero as (3 — > oo. For L > 6 the points fall very approximately onto a 
single curve; for both observables the finite-size corrections are very small. Our curve 
for the energy coincides with that of Shrock and Tsai ||15|| . The specific heat does not 



show any signal of transition points: it also decays smoothly to zero as (3 — > oo. 

In conclusion, there is no signal of phase transitions at any temperature T > in 
the 3-state hexagonal-lattice Potts antiferromagnet. This model is disordered at all 
temperatures. 

From the energy density one can easily compute the entropy density (per spin) of 



this model 33 



S(f3) = 5(0) + (3E((3) - f E{(3') dp , (3.3) 



o 



7 "Confidence level" is the probability that \ 2 would exceed the observed value, assuming that 
the underlying statistical model is correct. An unusually low confidence level (e.g., less than 5%) 
thus suggests that the underlying statistical model is incorrect — the most likely cause of which 
would be corrections to scaling. 



8 



where the value of the entropy at (3 = is given by Sh c , q =3(P = 0) = logg = log 3. 

Using our numerical data, we are able to compute the value of the entropy at 
(3 = 9, and we have to extrapolate this value somehow to (3 = oo. One way to achieve 
this is the following: at very large (3 we are deep in the strong coupling limit, so we 
expect that the energy density behaves as E(/3) ~ Ae~^ + 0(e~ 2 ^).[| If this is the 
case, we can compute exactly the integral on the r.h.s. of fl3~3[) and relate the result 
to the energy density at (3: 

S(oo) = S(0) - E((3) - [ E({3') d{3' . (3.4) 

Jo 

Now the maximum value of (3 where we have computed the energy becomes a cutoff. 
For each lattice size L, we have computed the values of S(j3 = oo, L) with this method 
using different values of the cutoff (3; the results were consistent within statistical 
errors (usually the differences were much smaller than the statistical errors). The 
values of the entropy density at zero temperature are displayed in Table [l|. The error 
bars are the sum of the statistical errors (coming from the statistical errors of the 
energies) and the systematic errors coming from the integration algorithm. We have 
used several extended trapezoidal rules and different sizes of the integration intervals 



34 1 : the systematic error takes account (conservatively) of the differences we found. 
If we fit the data to a constant [= 5^3(00, 00)] we find a good fit only for L m i n = 24 

^,3(00, 00) = 0.506844 ± 0.000012 (3.5) 

with x 2 — 0.11 (1 DF, level = 74%). This number is in agreement with Shrock and 

S hcS {oo, oo) ST = 0.5068 ± 0.0003 , (3.6) 



Tsai 15, 35 



but the error bar is one order of magnitude smaller than in Ref. [35]. If we use the 



extended Ansatz of Ref. 35 



ShcAP = 00, L) = S hc , 3 (oo, 00) + jp, + JZ + || 7 ( 3 - 7 ) 
we can reasonably fit all data (i.e., L min = 3) giving 

S hc>3 {oo, 00) = 0.506841 ± 0.000018 , (3.8) 

with x 2 = 0-40 (2 DF, level = 82%). This estimator agrees within errors with 

D/fD-0 



8 Although there is no obvious way to perform low-temperature series expansion for this model 
- there are too many inequivalent ground states-, it is reasonable to expect that there is an expansion 
in powers of , which corresponds to the minimum energy cost for a "overturned" spin. Indeed, 
our numerical data behaves in this way for large enough [3. 

9 An Ansatz of the form Shc,d,{P = 00, L) = 5^.3(00,00) + ci/L 2 is unable to fit well the 
data for any value of L m i n . If we add a term C2/L 4 , we get a good fit for L m i n = 6, giving 



5 ftc , 3 (oo,oo) = 0.506847 ± 0.000013 with X 2 = 0.075 (1 DF, level = 78%). The Ansatz Q is 
the first one to be able to fit all the data (L m i n — 3). Adding a term C4/L 8 does not modify the 
conclusions. 
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Indeed, the estimator ( j3.8|) is contained within the rigorous lower and upper 
bounds obtained by Shrock and Tsai for the hexagonal lattice [^, |37|: 



(l^l+m^m±Wl < e W~,~) < _ 3 , + 3)V» . (3.9) 
q - 1 

For g = 3 these bounds become 

0.505800 . . . = log (^~) < S hCt3 (oo, oo) < \og(V3) = 0.549306 .... (3.10) 

The lower bound is remarkably sharp: in Ref. [37| it is shown that if we extract the 
leading term [= q(l — 1/q) 3 ^ 2 } and expand the rest in powers of y = l/(q — 1), the 
resulting series for e Shc '^ 00 ' 00 ' } and its rigorous lower bound [cf. ( [3.9|) 1 agree up to 
0(y 10 ). The lower bound gives a very good approximation even for q as low as q = 3. 

The zero-temperature entropy density (|3.8| ) is a large fraction 46%) of the 
entropy at T = oo [Shc,3(P = 0) = log3 = 1.09861...]. This large ground-state 
entropy makes the system disordered at zero temperature. 



Acknowledgments 

We wish to thank Alan Sokal for encouragement and illuminating discussions, 
Robert Shrock for correspondence, and Chris Henley for making available some of his 
unpublished notes. The author's research was supported in part by CICyT grants 
PB95-0797 and AEN97-1680. 



10 



References 

[1] R.B. Potts, Proc. Camb. Philos. Soc. 48, 106 (1952). 

[2] F.Y. Wu, Rev. Mod. Phys. 54, 235 (1982); 55, 315 (E) (1983). 

[3] F.Y. Wu, J. Appl. Phys. 55, 2421 (1984). 

[4] R. Kotecky, cited in H.-O. Georgii, Gibbs Measures and Phase Transitions (de 
Gruyter, Berlin-New York, 1988), pp. 148-149, 457. 

[5] J. Salas and A.D. Sokal, J. Stat. Phys. 86, 551 (1997). 

[6] R.J. Baxter, Exactly Solved Models in Statistical Mechanics (Academic Press, 
London-New York, 1982). 

[7] R.J. Baxter, H.N.V. Temperley and S.E. Ashley, Proc. R. Soc. London A 358, 
535 (1978). 

[8] R.J. Baxter, J. Phys. A 19, 2821 (1986). 

[9] R.J. Baxter, J. Phys. A 20, 5241 (1987). 

[10] H. Saleur, Commun. Math. Phys. 132, 657 (1990). 

[11] PW. Kasteleyn and CM. Fortuin, J. Phys. Soc. Japan 26 (Suppl.), 11 (1969). 

[12] CM. Fortuin and PW. Kasteleyn, Physica 57, 536 (1972). 

[13] CM. Fortuin, Physica 58, 393 (1972); 59, 545 (1972). 

[14] R.J. Baxter, S.B. Kelland and F.Y. Wu, J. Phys. A9, 397 (1976). 

[15] R. Shrock and S.-H. Tsai, J. Phys. A 30, 495 (1997). 

[16] N. Metropolis, A.W. Rosenbluth, M.N. Rosenbluth, A.H. Teller and E. Teller, 
J. Chem. Phys. 21, 1087 (1953). 

[17] R.H. Swendsen and J.-S. Wang, Phys. Rev. Lett. 58, 86 (1987). 

[18] J.-S. Wang and R.H. Swendsen, Physica 167A, 565 (1990). 

[19] R.J. Baxter, Proc. R. Soc. Lond. A 383, 43 (1982). 

[20] J.-S. Wang, R.H. Swendsen and R. Kotecky, Phys. Rev. Lett. 63, 109 (1989). 

[21] J.-S. Wang, R.H. Swendsen and R. Kotecky, Phys. Rev. B 42, 2465 (1990). 

[22] S.J. Ferreira and A.D. Sokal, Phys. Rev. B 51, 6727 (1995). 

[23] J. Salas and A.D. Sokal, "The 3-state Square-Lattice Potts Antiferromagnet at 
Zero Temperature" , |cond-mat /980 1 079| . 



11 



[24] S.J. Ferreira and A.D. Sokal, in preparation. 

[25] S. Caracciolo, A. Pelissetto and A.D. Sokal, Phys. Rev. Lett. 72, 179 (1994). 

[26] R.J. Baxter, J. Math. Phys. 11, 784 (1970). 

[27] J. Kondev and C.L. Henley, Nucl. Phys. B 52, 6628 (1996). 

[28] C.L. Henley, in preparation. 

[29] J.K. Burton Jr. and C.L. Henley, J. Phys. A 30, 8385 (1998). 

[30] N.W. Ashcroft and N.D. Mermin, Solid State Physics (Sounders, Philadelphia, 
1976). 

[31] M. Campostrini, A. Pelissetto, P. Rossi and E. Vicari, "The two-point correla- 
tion function of three-dimensional O(N) models: critical limit and anisotropy", 
|cond-mat/9705086 . 

[32] N. Madras and A.D. Sokal, J. Stat. Phys. 50, 679 (1988). 

[33] K. Binder, Zeit. f. Physik B 45, 61 (1981). 

[34] W.H. Press, S.A. Teukolsky, W.T. Vetterling and B.P. Flannery, Numerical 
Recipies in C: The Art of Scientific Computing (Cambridge University Press, 
1992). 

[35] R. Shrock and S.-H. Tsai, Phys. Rev. E 55, 5165 (1997). 
[36] R. Shrock and S.-H. Tsai, Phys. Rev. E 55, 6791 (1997). 
[37] R. Shrock and S.-H. Tsai, Phys. Rev. E 56, 2733 (1997). 



12 



L 


a ( t \ 
Ohc,3{00,L) 


Q 
O 


O.oooooi ± U.UUUloo 


6 


0.509207 ±0.000097 


9 


0.507220 ±0.000089 


12 


0.506952 ± 0.000043 


24 


0.506864 ±0.000062 


48 


0.506843 ±0.000012 



Table 1: Values of the 3-state hexagonal-lattice Potts antiferromagnet entropy density 
at zero temperature Sh c ,3{(3 = oo, L) as a function of the lattice size L. The error 
bars are the sum of the statistical and systematic errors. 
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Figure 1: Curves where the hexagonal-lattice Potts model has been solved: (e J — 
l) 3 — 3q(e J — 1) = q 2 (solid curve), which has three branches; and the line e J = 
1 — q (dashed line). The horizontal dotted lines correspond to e J = 1 (separating 
the ferromagnetic and antiferromagnetic regimes) and to e J = (separating the 
antiferromagnetic regime from the unphysical region e J < 0). The squares (□) show 
the known ferromagnetic critical points (q = 1,2,3,4); and the diamond (O) marks 
the known antiferromagnetic critical point for q = 2. The open circles (o) show 
the points where the two antiferromagnetic branches cross the e J = line, namely 
q=(3±V5)/2. 
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q=3 hexagonal — lattice AF Potts model 




Figure 2: Second- mo merit correlation length £ of the 3-state hexagonal-lattice Potts 
antiferromagnet as a function of (3 and L. Symbols indicate L = 3 (+), L = 6 (O), 
L = 9 (□), L — 12 (o), and L = 24 (x). Points with L = 48 coincide with L = 24, 
and they are not shown for clarity. The isolated points on the far right of the picture 
(displayed for convenience at j3 — 9.5) are the data from the runs at (5 — oo. 
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q=3 hexagonal — lattice AF Potts model 
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Figure 3: Staggered susceptibility Xstagg of the 3-state hexagonal-lattice Potts anti- 
ferromagnet as a function of (3 and L. Symbols are as in Figure |2|. 
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